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Abstract 

The calculation of the two-loop corrections to the three jet production rate and to event shapes 
in electron-positron annihilation requires the computation of a number of up to now unknown two- 
loop four-point master integrals with one off-shell and three on-shell legs. In this paper, we compute 
those master integrals which correspond to planar topologies by solving differential equations in the 
external invariants which are fulfilled by the master integrals. We obtain the master integrals as 
expansions in e = (4 — d)/2, where d is the space-time dimension. The results are expressed in terms 
of newly introduced two-dimensional harmonic polylogarithms, whose properties are shortly discussed. 
For all two-dimensional harmonic polylogarithms appearing in the divergent parts of the integrals, 
expressions in terms of Nielsen's polylogarithms are given. The analytic continuation of our results to 
other kinematical situations is outlined. 



1 Introduction 



The calculation of perturbative next-to-next-to-leading order corrections to 2 — > 2 scattering or 1 — > 3 decay 
processes is a yet outstanding task for many precision applications in particle physics phenomenology. 
Progress in this field was up to very recently hampered by difficulties in the calculation of the virtual 
two-loop corrections to the corresponding Feynman amplitudes. 

Using dimensional regularisation with d = 4 — 2e dimensions as regulator for ultraviolet and 

infrared divergences, the large number of different integrals appearing in the two-loop Feynman amplitudes 
for 2-^2 scattering or 1 — > 3 decay processes can be reduced to a small number of master integrals. The 
techniques used in these reduction are integration-by-parts identities |^-^] and Lorentz invariance H. A 
computer algorithm for the automatic reduction of all two- loop four-point integrals was described in||. 

For two-loop four-point functions with massless internal propagators and all legs on-shcll, which are re- 
levant for example in the ncxt-to-ncxt-to-leading order calculation of two jet production at hadron colliders, 
all master integrals have been calculated over the past year [f7|-|l2|] . To compute the next-to-next-to-leading 
order corrections to observables such as the three jet production rate in electron-positron annihilation, two 
plus one jet production in deep inelastic electron-proton scattering or vector boson plus jet production 
at hadron colliders, one requires a different class of integrals: two-loop four-point functions with massless 
internal propagators and one external leg off-shell. Since these functions involve one scale more than their 
on-shell counterparts, one expects the corresponding master integrals to be more complicated, and also 
to be more numerous. First progress towards the computation of these master integrals has been made 
very recently by Smirnov, who computed the planar double box integral with all propagators raised to unit 
power [fl3"|| . Using a Mellin-Barnes contour integral technique, Smirnov obtained an analytic expression for 
the divergent parts and a one-dimensional integral representation for the finite term. Complementary work 
on the purely numerical evaluation of this type of master integrals has been presented recently by Binoth 
and Heinrich fl4|| . 

In this paper, we compute all master integrals appearing in the reduction of planar two-loop four- 
point functions with one external leg off-shell. The method employed here is substantially different from 
the techniques used in [l4|]: for all master integrals under consideration, we derived inhomogeneous 
differential equations |(| in the external invariants. The master integrals are then determined by solving 
these equations in terms of newly introduced two-dimensional harmonic polylogarithms and subsequent 
matching of the solution to the boundary condition. 

In Section ||, we briefly review the use of the differential equation approach to the computation of master 
integrals. To solve the differential equations for two-loop four-point master integrals with one off-shell leg, 
we employ an ansatz which is the product of a rational function of the kinematic invariants times the sum 
of two-dimensional harmonic polylogarithms (2dHPL). We explain how the unknown parameters in the 
ansatz are determined from the differential equation and its boundary condition. To illustrate how the 
method works in practice, we outline the calculation of two master integrals in Section || 

The master integrals computed using this method are listed in Section ^, which also contains a reference 
list of simpler master integrals for planar two- and there-point functions, which appear in the reduction of 
the four-point functions considered here. Section [| contains concluding remarks and an outlook on potential 
applications and open problems. 

We enclose an Appendix which summarises the properties of the harmonic polylogarithms (HPL) |lE( ] 
and extends the underlying formalism towards the newly introduced two-dimensional harmonic polylog- 
arithms (2dHPL). In the Appendix, we also tabulate relations between the HPL/2dHPL and the more 
commonly known Nielsen's generalised polylogarithms (l(],[l7j . We then give the expression of all (ordinary 
and two-dimensional) harmonic polylogarithms appearing in the divergent parts of the master integrals 
computed in this paper in terms of generalised polylogarithms of suitable arguments. As a consequence, 
the harmonic polylogarithms of weight 4 appearing in the finite parts of the master integrals can all be 
represented as one-dimensional integrals over ordinary Nielsen's polylogarithms up to weight three. Their 
precise numerical evaluation can be obtained by a simple extension of the techniques of ||1 7|| . 



f 



2 Differential Equations for Master Integrals 



The use of differential equations for the computation of master integrals was first suggested in [fL8| as 
means to relate integrals with massive internal propagators to their massless counterparts. The method 
was developed in detail in [fl9|| , where it was also extended towards differential equations in the external 
invariants. As a first application of this method, the two-loop sunrise diagram with arbitrary internal 
masses was studied in [^0 . We have developed the differential equation formalism for two-loop four-point 
functions with massless internal propagators, three external legs on-shell and one external leg off-shell in Jro. 
We derived an algorithm for the automatic reduction by means of computer algebra (using FORM |2l)| 
and Maple p2|) of any two-loop four-point integral to a small set of master integrals, for which we derived 
differential equations in the external invariants. For all master integrals with up to five internal propagators, 
we solved the differential equations in a closed form for arbitrary number of space-time dimensions d. In 
the present work, we elaborate on these results by deriving Laurent expansions around e = for all master 
integrals derived in Q as well as all remaining planar master integrals with six or seven propagators. 

Four-point functions depend on three linearly independent momenta: p\, P2 and p^. In our calculation, 
we take all these momenta on-shell (pf — 0), while the fourth momentum P123 = (pi + P2 + P3) is taken 
off-shell. The kinematics of the four-point function is then fully described by specifying the values of the 
three Lorentz invariants S12 = {pi +P2) 2 , S13 = {pi + P3) 2 and S23 = {p2 + Ps) 2 ■ 

Expressing the system of differential equations obtained in Q for any master integral in the variables 
S123 = S12 + S13 + S23, y — S13/S123 and z — S23/S123, we obtain a homogeneous equation in S123, and 
inhomogeneous equations in y and z. Since S123 is the only quantity carrying a mass dimension, the corre- 
sponding differential equation is nothing but the rescaling relation obtained by investigating the behaviour 
of the master integral under a rescaling of all external momenta by a constant factor. 

Some of the master integrals under consideration do not depend explicitly on all three invariants, but 
only on certain combinations of them, thus being one-scale or two-scale integrals. The one-scale integrals 
can not be determined using the differential equation method, since the only non-trivial differential equation 
fulfilled by them is the homogeneous rescaling relation. These integrals are relatively simple and can all be 
computed using Feynman parameters |^3lE3]. The two-scale integrals fulfil, besides the rescaling relation, 
one inhomogeneous differential equation, which can be employed for their computation. These integrals 
have all been computed in [^|,^5| ^7), where they are given in a closed form for arbitrary e. 

The solution of differential equations for two-loop four-point functions with one off-shell leg in terms 
of hypergeometric functions, was discussed in detail in ||. Although this method yields at first sight very 
compact results for arbitrary e, it is not very well suited for practical applications, where an expansion 
around e = is needed. The major problem with this approach is the appearance of generalised hyper- 
geometric functions which can be expressed only in terms of multidimensional integrals. The evaluation 
of these multiple integrals to the required order of the e-expansion is a problem comparable to the direct 
evaluation of the master integral from its Feynman parameter representation. To solve the differential 
equations for the three-scale integrals, we choose therefore an approach different from the one pursued 
in§. 

The procedure employed here is as follows. In the y and z differential equations for the master integral 
under consideration, the coefficient of the homogeneous term as well as the full inhomogeneous term (coef- 
ficients and subtopologies) are expanded as a series in e. For the master integral, we make an ansatz which 
has the form 

~^ fci(y, z; S123, e)Wi(y, z; e) , (2.1) 



where the prefactor 7ti(y, z; S123, e) is a rational function of y and 2, which is multiplied with an overall 
normalisation factor to account for the correct dimension in S123, while Tii(y, z;e) is a Laurent scries in e. 
The coefficients of its e-expansion are then written as the sum of 2dHPL up to a weight determined by the 
order of the series: 



H t (y,z;e)= f -Ye r ' 



T n(z)+^2 Y T n ^ j (z)H(m :j ;y) 

i=l m,j£Vj(z) 



(2.2) 
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where the H(fhj\y) are 2dHPL and T n (z), T n ^(z) are the as yet unknown coefficients. This form is 
motivated by the expectation that two-loop four-point functions can diverge up to 1/e 4 , as obtained by 
considering soft and collinear limits of internal propagators p8| ; p is an - a priori unknown - integer 
number, accounting for the fact that some of the master integrals might have a lower superficial degree 
of divergence. Vj(z) is the set of all possible indices for 2dHPL of weight j (j-dimensional vectors made 
from the components (0, 1,1 — z, z)). The T n (z) and T„ ^(2) are unknown coefficients, which can contain 
in turn ordinary one-dimensional HPL depending on z. With the ansatz ( [2.2] ), the leading l/e 4 -singularity 
corresponds to p — and n = 0. 

The rational functions TZi(y, z; S123, e) can be determined from the homogeneous part of the differential 
equations in y and z by inserting only the constant n = term of Tii(y, z\ e). 

It turns out that for all topologies which contain a single master integral, only one lZ(y, z; S123, e), and 



consequently also only one 7i(y,z;e), are present, such that the sum in (2.1) is not needed. For both 
topologies containing two master integrals, only one lZ(y, z; S123, e) and Tl{y,z;e) per master integral are 
sufficient, provided that the basis of first and second master integral is chosen appropriately, as explained 
later in this Section. 

Having determined the prefactor lZ(y, z; S123, e) for the master integral under consideration, we rewrite 
the y differential equation for the master integral into a y differential equation for H.(y,z;e). From the 
e-expansion of the inhomogeneous term, one can read off the value of p, required to match the order of the 
Laurent series of the master integral to the Laurent series of the inhomogeneous term. 

By construction of the 2dHPL, their derivative with respect to y is straightforward (see Appendix 



_d_ 

dy 



-r-H(m;y) = f(m;y) , 
dy 

Hfamj-ny) = f(m\y)H{rhj-r,y) , (2.3) 



where the f(m;y) are nothing but the y-dependent denominators present in the differential equation. 
Inserting the right hand side of ( |2.2| ) into the y differential equation for H(y, z; e) and differentiating 
according to the above rules, the differential equation becomes a purely algebraic equation. Grouping the 
different inverse powers of y, (1 — y), (1 — y — z) and (y + z) and using the linear independence of the 
H(fhj;y), this algebraic equation can be translated into a linear system of equations for the coefficients 
T nt rhj (z) ■ The resulting system of typically several hundreds of equations can be solved for the T n ^Az) 
using computer algebra (here we employ the same algorithm as used already in for solving large numbers 
of integration- by-parts and Lorentz-invariance identities). 

The T n (z), which do not multiply any y-dependent function in the ansatz and can therefore not be 
determined in the above way from the differential equation, correspond to the boundary condition. They 
are determined by exploiting the fact that all planar master integrals as well as their derivatives are analytic 
in the whole kincmatical plane with the exception of the two branch points at y = and z — 0. As a 
consequence, any factor (1 — y), (1 — y — z) or {y + z) appearing in the denominator of the homogeneous 
term of the differential equations for a master integral can be used to determine the boundary condition in 
y = 1, y = (1 — z) or y = —z respectively: multiplying the differential equation with this factor and taking 
the limit where the factor vanishes, one obtains a linear relation between the master integral in this special 
kinematical point and its subtopologies, no longer involving any derivatives. 

Some special care has to be taken in the case of topologies with two master integrals. Our choice of 
basis for the master integrals in these cases is as follows. For the five-propagator (diagonal box) master 
integrals, we use the basis already employed in H: 



P123 




d d k d d l 1 



P123- 
Pi- 



-Pi 

-P3 



(2ir) d {2-kY k 2 (k- Pl23 ) 2 (k - l- P 2) 2 l 2 (l-P3)' 2 

d d k d d l 1 
(2w) d (2n) d k 2 (k- Pl23 ) 2 (k - l-p 2 yp{l- P 3f 



(2.4) 



(2.5) 
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For the seven-propagator double box integrals, we use: 



P123- 

Pi- 

Pl23- 

Px- 



-P2 



-P3 



d d k d d l 



1 



(2) 



■*~P2 
+-P3 



(2ir)d (27r)" fc 2 (fc -p 23 ) 2 (fc -Pi 23 ) 2 (fc - imi-p 2 ) 2 (l-p 23 y 
d d k d d l {k-p 2 ) 2 

(2ir) d (2ir) d k?(k - p 23 ) 2 (k - P i2 3 ) 2 (k - - p 2 )*(l - p 23 ) 2 



(2.6) 



(2.7) 



This choice of basis for the double box integrals has been first suggested by Anastasiou, Tausk and Tejeda- 
Yeomans in [[f2| for the double box integrals with all legs on-shell. We find it to be convenient in the 
off-shell case as well. Apart from yielding a compact expression for the second master integral, this choice 
also circumvents problems with subleading terms p{| appearing in the reduction of tensor integrals. 

For both topologies, one obtains coupled sets of differential equations for the two master integrals. If 
only the n = term of ?ii(y,z;e) is retained, the equations decouple, thus enabling the determination of 
the 1Zi(y, z; S123, e) for each of the two master integrals. The boundary conditions for both master integrals 
are determined as above by using analyticity properties in the kinematical variables. Boundary conditions 
in two kinematical points are needed to fully constrain both master integrals. In the case of the seven 
propagator double box integrals (2.6,|]7), only one of those can be obtained in y = (1 — z), no other 
y-dependent denominator corresponding to an analytic point is present in the homogeneous term. To find 
the second boundary condition, wc first determine the z-dependence of the T n (z) by investigating their z 
differential equations, obtained from the y and z differential equations of the master integrals in y = (1 — z). 
The remaining constant terms in T n {z) are then obtained by using the fact that the master integrals are 
regular in z = 1. We would like to point out that this procedure relies crucially on the possibility to vary 
y and z independently, which is not the case for the on-shell double box integral, where the kinematical 
properties of the master integrals are insufficient [pi]] to determine both boundary conditions from the 
differential equation alone. 

Following the steps described in this section, we have derived the e-expansions of all planar four-point 
master integrals with one off-shell leg. It is obvious that this calculation had to follow a bottom-up 
approach starting from master integrals with few different denominators, which subsequently appeared in 
the inhomogeneous term of the differential equations for master integrals with more different denominators. 

As a check on all three-scale integrals obtained with this method, the results were inserted into the 
corresponding z differential equation. After transforming the H(rhj G Vj(z);y) into H(rhj € Vj{y); z), the 
z derivatives can be carried out straightforwardly, thus providing a check on our results. 

In the following Section, we illustrate on the example of the master integrals (2.4,2.5) how our method 
is applied in practice. In Section ||, we list our results for all master integrals appearing in the reduction of 
any planar two-loop four-point integral with one off-shell leg. 



3 Computing a Master Integral from its Differential Equation 



The diagonal box master integrals ( p.4|j2.5| ) are symmetric under the interchange of p± <-> p 2 , and conse- 
quently symmetric under y «-> z. The coupled set of y differential equations for them reads: 



d_ 
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-P3 
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+ (d-4) 
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4 1 - y — z 
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3(d-4) 2 y + z 
4 yz(l-y-z) 

3(d-3)(3d- 10) -l + y + 2z 



Pi- 



-Pa 
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-P23 



4 yz{l-y-z) 
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2(d-4) y 2 (l-y)(l-y-z) 

3(d-3)(3d-8)(3d-10) 1 



-Pi 



Pl3 




2(d-4) 



yz 2 (l - y - z) 




(3.1) 



(3.2) 



From the power s of (d — 4) multiplying the inhomogeneous terms, it can be read off that the second 
master integral ( |2.5| ) is suppressed by one power of (d — 4) with respect to the first master integral ([Q|). 
Expanding the homogeneous term of the differential equations and retaining only the leading contribution, 
the homogeneous differential equations decouple: 



d_ 

dy 

d_ 

dy 



P123- 

Pl- 
Pl23- 
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-P2 
-P3 

-Pa 

-P3 
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P123- 



Pl- 



-Pa 



-P3 



Pl23- 



Pl- 



-P2 



(3-3) 



(3.4) 



-P3 



The homogeneous z differential equations follow from the above by y <-* z, using the y «-» z inter- 
change symmetry of the integrals. As a result, one finds the prefactors for the first master integral (^J): 
TZ{y, z; S123, e) ~ l/(y + z) and for the second master i ntegral (|1|): TZ(y, z; S123, e) ~ l/(yz). 

Separating off these prefactors, we can transform (34_ 3^) into a set of coupled dif ferential equations 
for the two H(y, z; e), corresponding to both master integrals. Inserting the ansatz (|2.2j ) for the H(y, z; e), 
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we can then de term ine by comparison with the e-expansion of the inhomogeneous term that p = for (2.4) 
and p = 1 for (2.5). The coefficients T n ^,{z) appearing in the ansatz are then determined as outlined in 
the previous Section. 

To dete rmine t he boundary terms T n (z) for both master integrals, we first investigate the differential 
equations (3.1,3.2) in the point y = 1 — z, where both master integrals are analytic. Multiplying both 
equations with (1 — y — z) and taking y = 1 — z, the left hand sides of the equations vanish. From (3.1) we 
obtain, 



= — 



d-4 



P123 




z(l-z) 



P123- 



Pl- 



-P2 



-P3 



(Subtopologies) , 



(3-5) 

while ([T^) yields a multiple of this. The point y = (1 — z) can therefore be used only to determine a linear 
combination of boundary terms from both master integrals. A second boundary condition is obtained by 
multiplying ([O]) with (y + z) and subsequently taking y = — z, which fixes the boundary terms for the 
first master integral (|2.4|). The boundary terms for the second master integral ( |2.5| ) follow then from (B.J:), 
thus completing the computation of both master integrals. The results are given in the following section, 
together with the results for the other master integrals. 



4 Master Integrals 



In this Section, we tabulate all master integrals relevant for the computation of planar two-loop four point 
functions with one off-shell leg. We classify the integrals according to the number of different kinematical 
scales on which they depend into one-scale integrals, two-scale integrals and three-scale integrals. 
The common normalisation factor of all master integrals is 



(4tt) 



: r(f + e )r 2 (f- £ ) 

T(l - 2e) 



(4.1) 



4.1 One-scale Integrals 

The one-scale two-loop integrals with massless internal propagators were computed a long time ago in 
the context of the two- loop QCD corrections to the photon-quark-antiquark vertex [^3|Q. We list these 
integrals here only for completeness: 
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4.2 Two-scale Integrals 

The relevant two-scale integrals, corresponding to three-point functions with two off-shell legs, fulfil one 
inhomogeneous differential equation in the ratio of the two external invariants. Following the strategy 
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outlined in Section we computed these master integrals by employing a product ansatz of the form (2.1), 
with 77 containing only one-dimensional harmonic polylogarithms. 

Results for all two-scale integrals existed already in the literature in a closed form for arbitrary e in 
terms of ordinary hypergeometric functions |6,25-27f|. We checked that our results, obtained from solving 
the differential equations for these integrals, agree with the e-expansions of the hypergeometric functions 
quoted in §,|§-§7|]. 

Our results for the two-scale integrals read: 
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(4.17) 

(4.18) 
(4.19) 



4.3 Three-scale Integrals 

It was outlined in Section [2] that each three-scale master integral fulfils two inhomogeneous differential 
equations. For the computation of the master integral, it is sufficient to solve one of the differential 
equations, the second one is then used as an independent check on the result. Following the steps discussed 
in Section |^ we have computed all three-scale master integrals corresponding to planar topologies. 

The two different bubble insertions into the one-loop box integral with one off-shell leg were first 
considered in (27), where they were computed using the negative dimension approach. We confirmed 
these results in the differential equation approach. Both |2^| and || quoted a result in terms of generalised 
hypergeometric functions, valid for arbitrary e. As explained above, there is no straightforward procedure 
to expand these generalised hypergeometric functions in a Laurent series in e. Performing the Laurent series 
in e on the differential equations yields the e-expansion of the master integrals without the need to expand 
generalised hypergeometric functions. We obtain for bubble insertions into the one-loop box integral: 



P123- 

Pi- 



-Pi 
-P3 



16tH 



2 / N-2e 3 f K 1 . ( -SIS- 

{~ s 123) ' V Sl23 ' Sl2a 



S12 + S13 



0(e 2 ), 



(4.20) 



with: 



h.i,z{y,z) 
h.x,2{y,z) 

h.i,i{y,z) 

h.\,o{V' z ) 



= 



TV" 

~6 



-H(0;z) , 

-H(0; y)H(0; z) - 2iJ(0; z) + 2H(0, 0; z) + H(l, 0; z) 

-2H(0; y)H(0; z) - 2H(0; y)H(l, 0; z) - 4# (0; z) - H(0: z)H{\ - z, 0; y) 

-2H(Q, 0; y)H(0; z) + 4ff(0, 0; z) - 2H(0, 0; z)H(0; y) - 4F(0, 0, 0; z) - H(0, 1, 0; y) 

-2H(0, 1, 0; z) + 2H(1, 0; z) - H(l, 0; z)H(l - z;y) - 2H(1, 0, 0; z) - 2H(1, 1, 0; z) 



(4.21) 
(4.22) 



(4.23) 



H{1 - z, 1, 0; y) + — [+2 - 2H(0; y) - 3H(0; z) - 2H(l;z) - H(l - z; y)] 
b 



(4.24) 



/ 5 .i,-i(y, z) = +4H(0; y)H(0; z) - 4H (0; y)H(l, 0; z) + 4iJ(0; y)H(l, 0, 0; z) + 4H(0; y)H(l, 1, 0; z) 
-8iJ(0; z) - 2H(0; z)H{\ - z, 0; y) + 2H(0; z)H{\ - z, 0, 0; y) 
+H(0; z)H(l - z, 1 - z, 0; y) - 4ff(0, 0; y)H(0; z) + AH(0, 0; y)H(0, 0; z) 
+4H(0, 0; y)H(l,0; z) + 8H(0, 0; z) - 4ff(0, 0; z)H(Q; y) + 2H(0, 0; z)H(l - z, 0; y) 
+4ff(0, 0, 0; y)H(0; z) - 8H(0, 0, 0; z) + 4ff(0, 0, 0; z)H(0; y) + 8H(0, 0, 0, 0; z) 
+2H(0, 0, 1, 0; y) + 4F(0, 0, 1, 0; z) - 2F(0, 1, 0; y) - 4#(0, 1, 0; z) 
+4ff(0, 1, 0; z)ff(0; y) + 27J (0, 1, 0; z)#(l - z; y) + 2H(Q, 1, 0, 0; y) + 4ff(0, 1, 0, 0; 2) 
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-H(0, 1, 1, 0; y) + 4if (0, 1, 1, 0; z) + 2H(0, 1 - z; y)H(l, 0; z) 

+2if (0, 1 - z, 0; y)if (0; z) + 2if (0, 1 - z, 1, 0; y) + 4if (1, 0; z) - 2if (1, 0; z)H{\ - z; y) 
+27/(1, 0; z)H{\ - z, 0; y) + if(l, 0; z)if (1 - z, 1 - z; y) - 4if (1, 0, 0; z) 
+2if(l,0,0;z)if(l - z;y) + 4^(1,0,0,0; z) + 4ff(l, 0,1,0; z) - 4if (1, 1,0; z) 
+2ff(l, 1, 0; z)if(l - z; y) + 4if (1, 1, 0, 0; z) + AH {1, 1, 1, 0; z) + if (1 - z, 0, 1, 0; y) 
-2if (1 - z, 1, 0; y) + 2if (1 - z, 1, 0, 0; y) - if (1 - z, 1, 1, 0; y) + if (1 - z, 1 - z, 1, 0; y) 

+^ + 5C 3J ff(0;z) 
7T 2 r 

+ _ + 4 - 4if (0; y) + 4if (0; y)if (0; z) + 4if (0; y)H(l;z) - 6if (0; z) 
6 L 

+if (0; z)H(l - z; y) + 4if (0, 0; y) + 6if (0, 0; z) - if (0, 1; y) + 4if (0, 1; z) 
+21/(0, 1 - z; y) - 4if (1; z) + 2if (1; z)H{\ - z; y) + 4if (1, 0; z) + 4if (1, 1; z) 

-2if (1 - z; y) + 2JT(1 - z, 0; y) - if (1 - z, 1; y) + if (1 - z, 1 - z; y)l . (4.25) 



Pl23-»»- 
f>l-«- 



/ c 1 \2 / a N-2e 3 f 52 JjLUL -SZL.) 



\167T 2 / S23 
-P3 *=- 1 



with: 

/b.2, 3 (|/,*) = -1, (4-27) 
h.2,i{y,z) = -2 + H(0:y) + H(0;z) , (4.28) 
/ 5 .2,i(j/, z) = -4 + 2JT(0; y) - if (0; y)if (0; z) + 2if (0; z) - 2if (0, 0; y) - if (0, 0; z) - if (1, 0; y) , (4.29) 
/b.2,o(», «) = -8 + 4if (0; y) - 2if (0; y)/f (0; z) + if (0; y)if (1, 0; z) + 4if (0; z) + if (0; z)if (1 - z, 0; y) 
-4if (0, 0; y) + 2if (0, 0; y)H(0; z) - 2H(0, 0; z) + if (0, 0; z)H (0; y) + 4if (0, 0, 0; y) 
+if (0, 0, 0; z) + 2if (0, 1, 0; y) - 2if (1, 0; y) + if (1, 0; z)H(l - z; y) + 2if (1, 0, 0; y) 

+if (1, 1, 0; z) + if (1 - z, 1, 0; y) + 5( 3 + \ [+H (0; y) + if (1; z) + if (1 - z; y)] , (4.30) 

/5.2,-ifo, *) = -16 + 8if (0; y) - 4if (0; y)H(0; z) + 2if (0; y)H(l, 0; z) - if (0; y)if (1, 0, 0; z) 

-if (0; y)H(l, 1, 0; z) + 8if (0; z) + 2if (0; z)if (1 - z, 0; y) - 2if (0; z)H (1 - z, 0, 0; y) 

-if (0; z)if (1 - z, 1 - z, 0; y) - 8if (0, 0; y) + 4if (0, 0; y)H (0; z) - 2if (0, 0; y)if (0, 0; z) 

-2if (0, 0; y)if (1, 0; z) - 4if (0, 0; z) + 2if (0, 0; z)if(0; y) - if (0, 0; z)if (1 - z, 0; y) 

+8if (0, 0, 0; y) - 4if (0, 0, 0; y)H (0; z) + 2if (0, 0, 0; z) - if (0, 0, 0; z)H (0; y) 

-8if (0, 0, 0, 0; y) - if(0, 0, 0, 0; z) - 4if (0, 0, 1, 0; y) + 4if (0, 1, 0; y) 

-if (0, 1, 0; z)H (0; y) - if (0, 1, 0; z)if (1 - z; y) - 4if(0, 1, 0, 0; y) - if (0, 1, 1, 0; z) 

-2if (0, 1 - z; y)H(l, 0; z) - 2if (0, 1 - z, 0; y)if (0; z) - 2if (0, 1 - z, 1, 0; y) 

-4if (1, 0; y) + 2if (1, 0; z)H(l - z; y) - if (1, 0; z)H(l - z, 0; y) 

-if (1, 0; z)if (1 - z, 1 - z; y) + 4ff (1, 0, 0; y) - if (1, 0, 0; z)if (1 - z; y) 

-4if(l, 0, 0, 0; y) - ii(l, 0, 1, 0; z) + 2if(l, 1, 0; z) - ii(l, 1, 0, 0; z) 
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-2H(1 - z, 0, 1, 0; y) + 2H{\ - z, 1, 0; y) - 2H(1 - z, 1, 0, 0; y) 



-H(i- z ,i- z,i,0;y) + 



37tt 4 
360 



+ Ca [10 - 6H(0; y) - 5H(0; z) - H{1; z) - H(l - z; y)] 



+ 2H(0; y) - H(0; y)H(0; z) - H(0; y)H(l; z) - H(0; z)H{\ - z; y) 
2H(0, 0; y) - H(Q, 1; z) - 2H(0, l-z;y) + 2H(l; z) 
H(l, 0; z) + 2H{1 - z; y) - H{\ - z, 0; y) -H(l-z,l- z; y) . 



(4.31) 



The diagonal propagator insertion into the one-loop box with one off-shell leg gives rise to two different 
two-loop topologies, depending on whether the diagonal propagator is attached to the off-shell leg or not. 
One of these topologies contains one master integral, the other topology two master integrals. The choice 
of basis of master integrals for the latter and the separation of the differential equations are discussed in 
Sections || and ||[ 

The master integrals for both topologies were computed in terms of generalised hypergeometric functions 
valid for arbitrary e in Q . Again, an expansion in e of these functions is not straightforward, and we choose 
therefore to obtain the e expansions directly from the differential equations: 



P123- 



Pi- 



-Pi 



-P3 



16tt 2 



f \-2e 4 f.<>.( -213- _S23_ 

\— S123J ' V Sl23 ' Sl23 



S12 



+ 0(e), 



i=0 



with: 

/5.3,4(y,z) 
/5.3,3(y,2) 

f5.3,2(V' Z ) 

f5.3,i(y,z) 



, 
, 



(4.32) 



(4.33) 
(4.34) 



-H(0; y)H(0; z) - H(l, 0; y) - H(l, 0;z) - — , 

6 



(4.35) 

+2H(0; y)H(l, 0; z) + 2H(0; z)H(l - z, 0; y) + 2H(0, 0; y)H (0; z) + 2H(0, 0; z)H(0; y) 
+2H(Q, 1, 0; y) + 2H{0, 1, 0; z) + 2H(l, 0; z)H{l - z;y) + 2H(1, 0, 0; y) + 2H(1, 0, 0; z) 
+2H(1, 1, 0; z) + 2H(l - z, 1, 0; y) 



-— [+2H(0; y) + 2H(0; z) + 2H(1; z) + 2H(1 - z; y)} , 
6 



f5.3fl(Vi z ) 



(4.36) 

-4iJ(0; y)H(l, 0, 0; z) - AH(0; y)H(l, 1, 0; z) - 4iT(0; z)H{l - z, 0, 0; y) 

-4H(0; z)H(l - z ,l-z,0;y)- 4H(0, 0; y)H(0, 0; z) - 4ff (0, 0; y)if (1, 0; z) 

-4if (0, 0; z)ff (1 - z, 0; y) - 4#(0, 0, 0; y)ff(0; z) - 4#(0, 0, 0; z)iJ(0; y) 

-4iJ(0, 0, 1, 0; y) - AH (0, 0, 1, 0; z) - 4iJ(0, 1, 0; z)iJ(0; y) - AH(0, 1, 0; z)ff (1 - z; y) 

-4ff(0, 1, 0, 0; y) - AH(0, 1, 0, 0; z) - 4iJ (0, 1, 1, 0; z) - 4#(0, 1 - z; y)tf (1, 0; z) 

-4if (0, 1 - z, 0; y)tf (0; z) - 4if (0, 1 - z, 1, 0; y) - 41/(1, 0; z) if (1 - z, 0; y) 

-4ff (1, 0; z)H(l - z, 1 - z; y) - 4#(1, 0, 0; z)ff (1 - z; y) - 4iJ(l, 0, 0, 0; y) 

-AH(1, 0, 0, 0; z) - 4#(1, 0, 1, 0; z) - 4ff(l, 1, 0; z)H(l - z;y) - AH (1, 1, 0, 0; z) 

-4iJ(l, 1, 1, 0; z) - AH{1 - z, 0, 1, 0; y) - AH(1 - z, 1, 0, 0; y) - AH{1 - z, 1 - z, 1, 0; y) 

~ + y [ - 4#(0; y)iT(0; z) - 4H(0; y)ff(l; z) - 4#(0; (1 - z; y) - 4iT(0, 0; y) 
-47J(0, 0; z) - 4iJ(0, 1; z) - 4i? (0, 1 - z; y) - 4iJ(l; z)H(l - z;y) - 4ff(l, 0; z) 
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-4H(1, 1; z) - AH{1 - z, 0; y) - AH{1 -z,l-z;y) . 

P123-*- 



Pl- 



(4.37) 
(4.38) 



with: 

h.i,i(y,z) 
h.i,s{y,z) 

f5A,2(y,z) 

hA,i{y,z) 



hA,o{y,z) 



, (4.39) 

, (4.40) 

, (4.41) 

+H(0; y)H(l, 0; z) + H(0; z)H(l - z, 0; y) - H(0, 1, 0; y) - (0, 1, 0; z) 
+H(1, 0; z)ff (1 -z;y) + H(l, 1, 0; z) + (1 - z, 1, 0; y) 

+ ^[+tf(l;z) + tf(l-z;y)] , (4.42) 

-2H(0; y)H(l, 0, 0; z) - H(0; y)H(l, 1, 0; z) - 2tf (0; z)#(l - z, 0, 0; y) 
-if (0; z)H(l - z ,l-z,0;y) + 4H(0; z)H(z, 1 - z, 0; y) - 2JT(0, 0; y)tf (1, 0; z) 
-2if (0, 0; z)JT(l - z, 0; y) + 2JT(0, 0, 1, 0; y) - 2H(0, 0, 1, 0; z) + 2H(0, 1, 0; z)#(0; y) 
-2tf (0, 1,0; z)tf (1 - z; y) - 4tf (0, 1, 0; z)H(z; y) + 2H(0, 1, 0, 0; y) + 2£T(0, 1, 0, 0; z) 
-#(0, 1, 1, 0; y) + ff(0, 1, 1, 0; z) - 2tf (0, 1 - z; y)ff (1, 0; z) - 2H(0, 1 - z, 0; y)tf (0; z) 
-2H(0, 1 - z, 1, 0; y) - #(1, 0; z)H(l - z, 0; y) - H(l, 0; z)ff (1 - z, 1 - z; y) 
+4ff(l, 0; z)i?(z, 0; y) + 4if(l, 0; z)ff(z, 1 - z; y) - 2i? (1, 0, 0; z)ff(l - z; y) 
-2if(l, 0, 1, 0; z) + 4ff (1, 1, 0; z)H (z; y) - 2H(1, 1, 0, 0; z) - 2if(l - z, 0, 1, 0; y) 
-2if(l - z, 1, 0, 0; y) + if(l - z, 1, 1, 0; y) - if (1 - z, 1 - z, 1, 0; y) - AH{z, 0, 1, 0; y) 
+4H(z, 1 - z, 1, 0; y) + Cs [-2H (1; z) - 2H(1 - z; y)] 

7T 2 



- if (0; y)H(hz) - H(0; z)H{\ - z;y) - H(0, 1; y) + if(0, 1; z) - 2H(0, l-z;y) 
+4H(l; z)H(z; y) - H(l, 0; z) - if (1 - z, 0; y) + #(1 - z, 1; y) - if (1 -z,l-z;y) 



+4if(z, 1 - z; y) 



Pl23- 



Pi- 



-P2 



-P3 



= (J L _ 

V 16tt 2 



2 / \-2e 3 f, - . ( J13. 

(-S123) V"^ J5 - 5 ' 1 V«123 ' 



S23 
S123 



Sl3«23 



+ 0(e% 



i=-l 



with: 



fc.5,2(y>z) = -2H(0;y)-2H(0;z) 



+4if (0; y)if (0; z) + 4if (0, 0; y) + 4if (0, 0; z) + 3if (1, 0; y) + 3if (1, 0; z) + y 



(4.43) 



(4.44) 



(4.45) 
(4.46) 

(4.47) 



fa.5,i{y> z ) 

/ 5 .5,o(y, «) - -6H(0; y)H(l, 0; z) - 6if (0; z)H(l - z, 0; y) - 8ff (0, 0; y)if (0; z) - 8if (0, 0; z)/f (0; y) 
-8if (0, 0, 0; y) - 8if (0, 0, 0; z) - 6if (0, 1, 0; y) - 6if (0, 1, 0; z) - 6if (1, 0; z)if (1 - z; y) 
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-6/1(1, 0, 0; y) - 6/1(1, 0, 0; z) + 3/7(1, 1, 0; y) - 3/7(1, 1, 0; z) - 6/1(1 - z, 1, 0; y) 



h.t>-i{y,z) 



- 16 C3 + [-4H(0; y) - 4/1(0; z) + 31/(1; y) - 3//(l; z) - 67/(1 - z; y)] , 
= +12/7(0; y)//(l, 0, 0; z) + 6/7(0; y)H(l, 1, 0; z) + 12/7(0; z)/7(l - z, 0, 0; y) 
+67/(0; z)/7(l - z, 1 - z, 0; y) + 16/7(0, 0; y)/7(0, 0; z) + 12//(0, 0; y)H(l, 0; z) 
+12/7(0, 0; z)H{\ - z, 0; y) + 16/7(0, 0, 0; y)//(0; z) + 16/7(0, 0, 0; z)//(0; y) 
+16/7(0, 0, 0, 0; y) + 16/7(0, 0, 0, 0; z) + 12/7(0, 0, 1, 0; y) + 12/7(0, 0, 1, 0; z) 
+12/7(0, 1, 0; z)/7(0; y) + 12/7(0, 1, 0; z)H{\ - z; y) + 127/(0, 1, 0, 0; y) 
+12/7(0, 1, 0, 0; z) - 6H (0, 1, 1, 0; y) + 6/7(0, 1, 1, 0; z) + 121/(0, 1 - z; y)//(l, 0; z) 
+ 121/(0, 1 - z, 0; y)H(0; z) + 121/(0, 1 - z, 1, 0; y) + 61/(1, 0; z)H(l - z, 0; y) 
+611(1, 0; z)H{\ - z, 1 - z; y) + 12/1(1, 0, 0; z)//(l - z; y) + 12/7(1, 0, 0, 0; y) 
+ 12//(l,0,0,0;z) + 12//(1, 0,1,0; z) - 6//(l, 1, 0, 0; y) + 6//(l, 1, 0, 0; z) 
+3/1(1, 1, 1, 0; y) + 3/1(1, 1, 1, 0; z) + 12//(1 - z, 0, 1, 0; y) + 12//(1 - z, 1, 0, 0; y) 

2tt 4 

-6/1(1 - z, 1, 1, 0; y) + 6/1(1 - z, 1 - z, 1, 0; y) - — 

5 

+C 3 [+321/(0; y) + 32/1(0; z) + 31/(1; y) + 15/1(1; z) + 12/1(1 - z; y)] 



(4.48) 



7T 

IT 



8//(0; y)//(0; z) + 6/1(0; y)/Z(l; z) + 6Z/(0; z)/Z(l - z; y) + 8/Z(0, 0; y) 



+8/1(0, 0; z) - 6Z/(0, 1; y) + 6/1(0, 1; z) + 121/(0, 1 - z; y) + 6Z/(1, 0; z) + 3/1(1, 1; y) 



+31/(1, 1; z) + 6/1(1 - z, 0; y) - 61/(1 - z, 1; y) + 6Z/(1 - z, 1 - z; y) . 
The only six propagator planar master integral was up to now not known in the literature: 

1>\ 13 -^-| 1 

' E — v .v v -^). 



pi- 



-P3 



(JLV (- gl23 )- 2e ^ / 6 .m(S'S) 

\ 167T 2 / (si2 + s 23 )s23 4r! 



with: 

f6.iAy> z ) 
h.iAVi z ) 

fe.iAV'Z) 
h.\,i{y,z) 



= 



+/7(0,0;y)-/7(l,0;y) 



(4.49) 
(4.50) 

(4.51) 
(4.52) 

(4.53) 



fe.i,o(y,z) 



= +1/(0; y)H(l, 0; z) + 1/(0; z)//(l - z, 0; y) + Z/(0, 0; z)/7(0; y) - 21/(0, 0, 0; y) 

-1/(0, 1, 0; y) + //(0, 1, 0; z) + //(l, 0; z)H(l - z; y) + 21/(1, 0, 0; y) - 2//(l, 1, 0; y) 
+//(l,l,0;z) + //(l-z,l,0;y)-3C3 

+ ^ [+H(0: z) - 2//(l: y) + 1/(1; z) + H(l - z; y)] , (4.54) 

= -3/1(0; y)H(l, 0, 0; z) - 1/(0; y)//(l, 1, 0; z) + 21/(0; z)H(l, 1 - z, 0; y) 

-2/1(0; z)H(l - z, 0, 0; y) - 1/(0; z)/Z(l - z, 1 - z, 0; y) - 21/(0, 0; y)H(0, 0; z) 
-21/(0, 0: y)H(l, 0; z) + 21/(0, 0; z)//(l, 0; y) - 3/1(0, 0; z)H(l - z, 0; y) 
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-3H(0, 0, 0; z)H(0; y) + 4ff(0, 0, 0, 0; y) + 2H(0, 0, 1, 0; y) - 3H(0, 0, 1, 0; z) 

-if (Q, 1, 0; z)if (0; y) + 2H(0, 1, 0; z)if (1; y) - 3if (0, 1, 0; z)H(l - z;y) + 2if (0, 1, 0, 0; y) 

-3if (0, 1, 0, 0; z) + if (0, 1, 1, 0; y) - 2H(0, 1 - z; y)H(l, 0; z) - 2H(0, 1 - z, 0; y)H(0; z) 

-2if (0, 1 - z, 1, 0; y) + 2if (1, 0; y)H(l, 0; z) + 2if (1, 0; z)if (1, 1 - z; y) 

-if(l,0;z)if(l - z,0; y) - if(l, 0; z)if(l - z, 1 - z;y) - 3if (1, 0, 0; z)H (1 - z;y) 

-4if (1, 0, 0, 0; y) - 2if (1, 0, 1, 0; y) - 3if (1, 0, 1, 0; z) + 2H(1, 1, 0; z)if (1; y) 

+411 (1, 1, 0, 0; y) - 3H(1, 1, 0, 0; z) - 4H(1, 1, 1, 0; y) + 2H{1, 1 - z, 1, 0; y) 

-2H(1 - z, 0, 1, 0; y) - 2if(l - z, 1, 0, 0; y) + 2H(1 - z, 1, 1, 0; y) - H (1 - z, 1 - z, 1, 0; y) 



^ 4 




72 


+ 






6 





if(0; y)if (0; z) - If (0; y)H(l; z) + 2if (0; z)H(l; y) - if (0; z)H{\ - z; y) 
-H(0, 0; z) + if (0, 1; y) - 2H(0, 1 - z;y) + 2if (l;y)if(l; z) - #(1,0; z) - 4if(l, 1; y) 



-2fi(l, 1 - z; y) - H(l - z, 0; y) + 2H{1 - z, 1; y) - If (1 - z, 1 - z; y) 



(4.55) 



The seven propagator planar double box topology contains two master integrals. We choose the basis 
first suggested in [Q for the double box integrals with all legs on-shell, which yields particularly compact 
results also in our case with one leg off-shell: 



with: 

/7.i,2(y,z) 



Pl23- 
Pl- 



-Pl 
-P3 



I 16^ 2 / 



s 123 ) ' V S123 ' S123 



S13S| 3 



0(e), 



i=0 



2[if(0;y) + if(0;z)] , 

-4if (0; y)ff (0; z) - 4ff (0, 0; y) - 4ff(0, 0; z) - if (1, 0; y) - 3if (1, 0; z) - 

+6if (0; y)#(l, 0; z) + 2if (G; z)if (1 - z, 0; y) + 8if (0, 0; y)H(0; z) 
+8H(0, 0; z)if (0; y) + 8if (0, 0, 0; y) + 8if (0, 0, 0; z) + 5if (0, 1, 0; y) 
+6if(0, 1, 0; z) + 211(1, 0; z)if (1 - z;y) + 2if (1, 0, 0; y) + 9if (1, 0, 0; z) 

-2if (1, 1, 0; y) + 6JT(1, 1, 0; z) + 2if (1 - z, 1, 0; y) - % 



(4.56) 



(4.57) 
(4.58) 

(4.59) 



+ — [+7H(0; y) + 71/(0; z) - 21/(1; y) + 611(1; z) + 211(1 - z; y)] , 

A.i.oG/, z) = -1811(0; y)H(l, 0, 0; z) - 1211(0; y)H(l, 1, 0; z) + 211(0; z)H(l, 1 - z, 0; y) 

-411(0; z)H{\ - z, 0, 0; y) - 211(0; z)H(l - z, 1 - z, 0; y) - 1611(0, 0; y)H(0, 0; z) 
-1211(0, 0; y)H(l, 0; z) + 211(0, 0; z)H(l, 0; y) - 611(0, 0; z)H(l - z, 0; y) 
-1611(0, 0, 0; y)H(0; z) - 1611(0, 0, 0; z)if (0; y) - 1611(0, 0, 0, 0; y) 
-1611(0, 0, 0, 0; z) - 1611(0, 0, 1, 0; y) - 1211(0, 0, 1, 0; z) - 1211(0, 1, 0; z)H(0; y) 
+211(0, 1, 0; z)H(l; y) - 611(0, 1, 0; z)H(l - z; y) - 1011(0, 1, 0, 0; y) 



(4.60) 
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-1877(0, 1, 0, 0; z) + 77(0, 1, 1, 0; y) - 12/7(0, 1, 1, 0; z) - 1077(0, 1 - z; y)H(l, 0; z) 
-1077(0, 1 - z, 0; y)/7(0; z) - 1077(0, 1 - z, 1, 0; y) + 2H(1, 0; y)/7(l, 0; z) 
+2H(l, 0; 2)77(1, 1 - z; y) - 4/7(1, 0; z)H(l - z, 0; y) - 2H(1, 0; z)/7(l -z,l-z;y) 
-677(1, 0,0;z)77(l-z;y) -4/7(1, 0,0,0; y) - 21/7(1, 0, 0, 0; z) - 277 (1, 0, 1, 0; y) 
-1577(1, 0, 1, 0; z) + 2H(1, 1, 0; z)/7(l; y) - 477(1, 1, 0; z)77 (1 - z; y) 
+477(1, 1, 0, 0; y) - 2177(1, 1, 0, 0; z) - 4/7(1, 1, 1, 0; y) - 1277(1, 1, 1, 0; z) 
+2/7(1, 1 - z, 1, 0; y) - 2/7(1 - z, 0, 1, 0; y) - 4/7(1 - z, 1, 0, 0; y) 



+4/7(1 - z, 1, 1, 0; y) - 2/7(1 - z, 1 - z, 1, 0; y) - — 

+C 3 [+7/(0; y) - 2H(1 - z;y) - 6/7(1; y) + 3/7(1; z) + 77(0; z) 



TV 

~6 



- 1477(0; y)/7(0; z) - 12H (0; y)H(l; z) + 277(0; z)H(l; y) - 277(0; z)H(l - z; y) 

-1477(0, 0; y) - 14/7(0, 0; z) + 77(0, 1; y) - 1277(0, 1; z) - 1077(0, 1 - z; y) 
+277(1; y)77(l; z) - 477(1; z)77(l - z; y) - 1577(1, 0; z) - 477(1, 1; y) - 1277(1, 1; z) 



+277(1, 1 - z; y) - 477(1 - z, 0; y) + 477(1 - z, 1; y) - 277(1 - z, 1 - z; y) 



P123- 



(2) 



-P2 



S123J 



-2e 4 



_ p3 V167TV («ia + S13)S 23 l=0 



E " v T" 1M; +o( e ), 



(4.61) 
(4.62) 



= 0, 



O c: c_2 

+277(0; y)77(0; z) + -77(0, 0; z) + -77(1, 0; z) + — 



(4.63) 
(4.64) 



(4.65) 



= -477(0; y)77(l, 0; z) - 477(0; z)77(l - z, 0; y) - 477(0, 0; y)77(0; z) - 677(0, 0; z)77(0; y) 
~ 77(0, 0, 0; z) - 477(0, 1, 0; y) - y 77(0, 1, 0; z) - 477(1, 0; z)77(l - z; y) 

~ 77(1, 0, 0; z) - ^77(1, 1, 0; z) - 477(1 - z, 1, 0; y) 



7T 

~6~ 



-477(0; y) - y77(0; z) - ^77(1; z) - 477(1 - z; y) 



(4.66) 



+ 1277(0; y)77(l, 0, 0; z) + 877(0; y)77(l, 1, 0; z) + 877(0; z)77(l - z, 0, 0; y) 

+ 1077(0; z)77(l - z, 1 - z, 0; y) + 1277(0, 0; y)77(0, 0; z) + 877(0, 0; y)/7(l, 0; z) 

+ 1277(0, 0; z)77(l - z, 0; y) + 877(0, 0, 0; y)77(0; z) + 1477(0, 0, 0; z)77(0; y) 

+ y 77(0, 0, 0, 0; z) + 1077(0, 0, 1, 0; y) + y 77(0, 0, 1, 0; z) + 1277(0, 1, 0; z)77(0; y) 

27 

+ 1277(0, 1, 0; z)77(l - z; y) + 877(0, 1, 0, 0; y) + —H(0, 1, 0, 0; z) - 277(0, 1, 1, 0; y) 
+ y 77(0, 1, 1, 0; z) + 1077(0, 1 - z; y)77(l, 0; z) + 1077(0, 1 - z, 0; y)77(0; z) 
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+lCLff(0, 1 - z, 1,0; y) + 8H(1,0; z)H(l - z, 0; y) + 10ff(l, 0; z)if(l - z, 1 - z; y) 

35 25 
+12H (1, 0, 0; z)ff(l - z; y) + — (1, 0, 0, 0; z) + —H(l, 0, 1, 0; z) 

+8H(1, 1, 0; z)tf (1 - z; y) + y (1, 1, 0, 0; z) + yif(l, 1, 1, 0; z) + 10H(1 - z, 0, 1, 0; y) 

377T 4 

+8H(1 - z, 1, 0, 0; y) - 2ff (1 - z, 1, 1, 0; y) + 10£T(1 - z, 1 - z, 1, 0; y) + — 
+C 3 [+4IT(0; y) + 6H (1; z) + 4ff(l - z; y)] 

2 

' + 10ff (0; y)ff (0; z) + 8H(0; y)H(l\ z) + 10ff(0; z)F(l - z; y) + 8ff (0, 0; y) 



7T 



35 21 
+ y (0, 0; z) - 2H(0, 1; y) + yff(0, 1; z) + 101/(0, 1 - z; y) + 8ff (1; z)ff (1 - z; y) 

+ y (1, 0; z) + y (1, 1; z) + &ff(l - z, 0; y) - 2if (1 - z, 1; y) + 101/(1 - z, 1 - z; y) 



(4.67) 



The first of the three-scale double box integrals (4.56), was recently computed by Smirnov |L3| using 
a Mellin-Barnes representation. We confirm this result. Since a comparison of the result of |D| with our 
result is not straightforward, we outline below the procedure we have employed. 

In |l3| ] analytic forms containing Nielsen's polylogarithms are given for all e-divergent terms of ( 4.56| ); 



these terms can be compared directly. For the finite part of (4.56), Smirnov provides in [ |13| a result 
involving a one-dimensional integral plus other terms containing only Nielsen's polylogarithms of non- 
simple arguments. The integrand of the one-dimensional integral involves products of logarithms with 
dilogarithms. 

To compare with our result, we rewrote the finite part quoted in |l3[ in terms of 2dHPL. The procedure 
employed was as follows. In the limit y — > 0, we determined the z-dependent coefficients of In 1 y (i = 1 ... 4) 
as well as the finite z-dependent term. The remaining terms, which are regular in y — ► were determined 
from the y-derivative of the result. Taking this derivative, one obtains a combination of rational fractions 
1/y, 1/(1 — y), 1/(1 — y — z) times Nielsen's polylogarithms of level 3. These polylogarithms are differentiated 
with respect to y again, resulting in (1/y, 1/(1 — y), 1/(1 — y — z)) times Nielsen's polylogarithms of level 
2, which are translated into 2dHPL of weight 2. The two y derivatives are undone by integration (keeping 
correct account of boundary terms by verifying the limit y — > at each stage). As a result, we obtain 
2dHPL of weight 4, as well as products of 2dHPL with ordina ry H PL of z. After converting the overall 
normalisation factor, we find agreement of [|l3| with our result ( p~56l ). 



5 Conclusion 

Two-loop four-point functions with one off-shell leg are an important ingredient to the calculation of next- 
to-next-to-leading order corrections to three jet production and related observables in electron-positron 
annihilation. By exploiting integration-by-parts |^-^| and Lorentz invariance || identities, one can express 
the loop integrals appearing in these functions as a linear combination of a small set of master integrals. 
These master integrals are scalar functions of the external invariants, their determination was up to now a 
major obstacle to further progress in next-to-next-to-leading order calculations. 

In the present paper, we compute all master integrals appearing in the reduction of planar two- loop four- 
point functions with massless internal propagators and one off-shell leg. The method employed here relies on 
the fact that all master integrals fulfil inhomogeneous differential equations in their external invariants. For 
the master integrals under consideration, these differential equations were derived in We determined 
the master integrals by solving the corresponding differential equations, starting from a product ansatz 
involving a rational function of the external invariants and a sum of newly introduced two-dimensional 
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harmonic poly logarithms. The two-dimensional harmonic poly logarithms are an extension of the harmonic 
polylogarithms of |^5| , built a la carte to suit the needs of our calculation. We describe the underlying 
formalism in the Appendix and tabulate the relation of two-dimensional harmonic polylogarithms to the 
better known generalised polylogarithms of Nielsen |l(|. For the divergent parts of the master integrals 
considered here, we find analytic expressions in terms of generalised polylogarithms of weight up to 3, while 
the finite part involves also 2dHPL of weight 4 which can be expressed as a one-dimensional integral over 
a combination of generalised polylogarithms of weight 3 with suitable arguments. 

The application of our results to different kinematic situations, as encountered in the next-to-next- 
to-leading order corrections to vector boson plus jet production at hadron colliders or two plus one jet 
production in deep inelastic scattering requires the analytic continuation of the two-dimensional harmonic 
polylogarithms outside the region 0<y<l,0<z<l — y, where they are real. An algorithm to derive 
these analytic continuation formulae is outlined in the Appendix. 

A yet outstanding task is the computation of the master integrals appearing in the reduction of am- 
plitudes with crossed topologies. The derivation of these from the corresponding differential equations 
should follow the procedure outlined in Section ^. The ansatz does however contain more than one rational 
function IZi, as already observed in the on-shell case |||l(J. Moreover, the determination of the boundary 
conditions is more involved since the non-planar topologies have in general three branch points: in y = 0, 
z — and y = 1 — z. This work will be reported in a separate publication [pl|. 
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A Harmonic Polylogarithms 

Harmonic polylogarithms (HPL) were introduced in fl5|| as an extension of the generalised polylogarithms 
of Nielsen Jl6| . They are constructed in such a way that they form a closed, linearly independent set under 
a certain class of integrations. We observe that the class of allowed integrations on this set can be extended 
a la carte by enlarging the definition of harmonic polylogarithms in order to suit the needs of a particular 
calculation. We will make use of this feature by generalising the one-dimensional HPL of [jl5| to two- 
dimensional harmonic polylogarithms (2dHPL), which appear in the solution of the differential equations 
for the three-scale master integrals discussed in this paper. 

In the following subsection, we will briefly review the properties of the ordinary HPL, and tabulate 
expressions for HPL up to weight 4 in terms of generalised polylogarithms. In another subsection, we 
will then extend this framework to 2dHPL, which can in general be expressed up to weight 3 in terms 
of generalised polylogarithms of suitable non-simple arguments. We tabulate all 2dHPL up to weight 
3. The 2dHPL of weight 4 are expressed as one-dimensional integrals over combinations of generalised 
polylogarithms and therefore can be evaluated numerically in a straightforward manner. 

A.l One-dimensional Harmonic Polylogarithms 

The one-dimensional HPL H(rh w ;x) is described by a w-dimensional vector rh w of parameters and by its 
argument x. w is called the weight of H . We briefly recall the HPL formalism: 

1. Definition of the HPL at w = 1: 

H(l;x) = -ln(l-x), 

H(0;x) ee lnx, 

H{-l;x) ee ln(l + x) (A.l) 
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and the rational fractions in x 

/(i; *) 



1 



/(0; 



1-x 
1 

a; 



such that 



2. For to > 1: 



/(-!;*) = ~r~T ' ( A - 2 ) 
1 + x 

d 

—H{a;x) = f(a;x) with a = +1,0,-1. (A.3) 



H(0,...,0;x) = ^-r\n w x, (A.4) 

H(a,b;x) ee / dx7(o;ar')if(6;x') , (A.5) 
./o 

which results in 

^H(a, b- x) = /(a; :r)F(&; x) . (A.6) 

This last relation is the main tool for verifying identities among different HPL. Such identities can be 
verified by first checking a special point (typically x = 0) and subsequently checking the derivatives. 
If agreement in the derivatives is not obvious, this procedure can be repeated until one arrives at 
relations involving only HPL with w = 1. 

3. The HPL fulfil an algebra (see Section 3 of jll|), such that a product of two HPL (with weights w% 
and W2) of the same argument a: is a combination of HPL of argument x with weight w — wi + W2, 

H(a;x)H(b;x) = ^ H(c;x), (A.7) 

c=a<Sb 

where a W b represents all mergers of a and b in which the relative orders of the elements of a and b 
are preserved. For example at wi — 1,W2 — 3, this relation reads: 

H(mi; x)H(rri2, m^, x) = +H(mi, rn,2, ms,m4; x) +-ff(m2,?ni,m3, r/14; a;) 

+H(m,2, m 3 , mi,m&] x) + H(m,2, ma, m.4, mi; x) . (A. 8) 

4. The HPL fulfil the integration- by-parts identities (see also Section 3 of E]) 

H{m\, . . . , m g ; a;) = H(mr, x)H(rn2, • • • , a;) — H(rri2, mi; x)H{m^, . . . , m q ; x) 

+ ... + {-iy +1 H(m q , . . . ,mi; x) . (A.9) 

5. The HPL are linearly independent. 



The product identities ( |A.7| ) and the integration-by-parts identities ( |A.9j ) involve different polyloga- 
rithms of the same weight w, as well as products of harmonic poly logarithms of lower weight. They can 
thus be used to express all polylogarithms of weight w in terms of a so-called 'minimal set' of weight w 
plus products of HPL of low er weight. 

It can be seen from (A.5), that the HPL of parameters (+1,0, —1) form a closed set under the class of 
integrations 

^(h>ih<'Th) H ^- (A - 10) 
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In the context of the present calculation, it turns out that integrals involving denominators 1/(1 + x) do 
never occur, the HPL of parameters (0, 1) are therefore sufficient. 

The HPL of parameters (0, 1) can all be expressed in terms of logarithms and Nielsen's polylogarithms. 
Up to level m = 4, they read: 



H(0;x) 
H(l;x) 




In a; , 

-ln(l -x) , 


H(0,0;x) 

#(0,1; x) 
H(l,0;x) 


= 
= 


|m 2 x, 
Li 2 (x) , 

— In x ln(l — x) — Li 2 (x) , 


H(l,l;x) 


= 


^ln 2 (l-x), 


#(0,0, 0;x) 

#(0,0, l;x) 
#(0,1, 0;x) 
#(0,1, l;x) 




= 
= 


i In 3 x , 
o! 

Li 3 (x) , 

-2Li 3 (x) + lnxLi 2 (x) , 
Si, 2 (x) , 


#(1,0, 0;x) 




— i ln(l — x) In 2 x — lnxLi 2 (x) + Li 3 (x) , 


H(l,0,l;x) 




-2Si, 2 (x) -ln(l-x)Li 2 (x) , 


#(1,1, 0;x) 


= 


Si, 2 (x) +ln(l -x)Li 2 (x) + ilnxln 2 (l - x) , 


#(1,1, l;x) 


= 


-^ln 3 (l-x), 


#(0,0,0,0;x) 

#(0,0,0,1; x) 
#(0,0,1,0; x) 
#(0,0, 1,1; x) 


= 

= 


^ln 4 x, 
Li 4 (x) , 

lnxLi 3 (x) — 3Li 4 (x) , 
S 2 , 2 (x) , 


H(0,1, 0,0; x) 




i In 2 x Li 2 (x) — 2 In x Li 3 (x) + 3 Li 4 (x) , 
2 


#(0, 1,0,1; x) 




-2S 22 (x) + -[Li 2 (x)l 2 , 


#(0,1, 1,0; a;) 




lnx,S 1 , 2 (x)--[Li 2 (x)] 2 , 


#(0,l,l,l;x) 




Si, 3 (aO , 


#(l,0,0,0;x) 




— i In 3 xln(l — x) — i In 2 xLi 2 (x) + ln(x) Li 3 (x) — Li 4 (x) , 


#(1,0, 0,1; x) 




-i[Li 2 (x)] 2 -ln(l-x)Li 3 (x) , 


#(1,0, 1,0; x) 




21n(l - x) Li 3 (x) - lnxln(l - x) Li 2 (x) - 21nxSi ;2 (x) + i[Li 2 (x)] 2 + 2S 2 , 2 (x) 


#(1,0, 1,1; or) 




-ln(l-x)S li2 (x)-3Si, 3 (x) , 
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77(1, 1, 0, 0; a;) = — In 2 x ln 2 (l — x) — ln(l — x) \A^{x) + lnx ln(l — x) Li 2 (a:) + lnx Si. 2 (a;) — S 2;2 (a;) , 

77(1, 1, 0, 1; ar) = - ln 2 (l - x) IA 2 {x) + 2 ln(l - x) Si, 2 (ai) + 3 Si, 3 (a?) , 

77(1, 1, 1, 0; x) = -- In xln 3 (l - a;) - ^ ln 2 (l - ar) Li 2 (a;) - ln(l - x) Si, 2 (x) - Si,a(x) , 

77(1,1,1,1;^) - ~ln 4 (l-a:). (A.ll) 



A. 2 Two-dimensional Harmonic Poly logarithms 



The generalisation from one-dimensional to two-dimensional HPL starts from ( A.1C ), which defines the 
class of integrations under which the HPL form a closed set. By inspection of the various inhomogeneous 
terms of the y differential equations for the three-scale master integrals discussed in this paper, we find 
that, besides the denominators 1/y and 1/(1 — y) also 1/(1 — y — z) and l/(y + z) appear. It is therefore 
appropriate to introduce an extension of the HPL, which forms a closed set under the class of integrations 

v f 1 I 1 1 \ 

dW-1 7-1 ; >-TT- (A.12) 

Such an extension of the HPL formalism can be made by extending the set of fractions by 



1-y-z 
1 

y + z 



f(z;y) = (A.13) 



and correspondingly the set of HPL at weight 1 by 

77(1 -z;y) = -lnll- 



y 



H(z;y) = ln(^). (A.14) 



Allowing (z, 1 — z) as components of the vector rh w of parameters, ( |A.5| ) does then define the extended 
set of HPL, which we call two- dimen sional harmonic poly logarithms (2dHPL). They form a closed set 
under integrations of the form ( A.12] ), as required. This closure was achieved by construction, showing 



that the HPL formalism of |15[| can be extended to suit the needs of a particular calculation. We do 
not consider —1 among the components of rh w , since they are not needed in the present context. The 
2dHPL fulfil the same algebra as the HPL ( |A.7| ), and they are linearly independent. The feature of linear 
independence distinguishes the 2dHPL from Nielsen's polylogarithms of non-trivial two- variable arguments 
(in the literature, see for example |30), there are large numbers of relations between Nielsen's polylogarithms 
with arguments which are rational functions of two variables). Let us observe here that the use of linearly 
independent functions protects against the danger of writing "hidden zeroes" (complicated expressions 
whose actual value is identically zero) , making it trivial at the same time to carry out those simplifications 
whose occurrence is among the main features of any analytical calculation. Finally, the integration-by-parts 



relation (A. 9) remains valid also for 2dHPL. 

Two-dimensional harmonic polylogarithms can be expressed in terms of generalised polylogarithms only 
up to weight 3. At weight 4, only some special cases relate to generalised polylogarithms. The relations at 



weight 2 read: 

H(0,l-z;y) = Li 2 (-^ J , 
77(0, z;y) = -Li 2 



y 
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H(l,l-z;y) = iln 2 (l-y)-ln(l-y)ln(l-z)+ Li 2 (j^j ~ Li 2 (z) 
H(l,z;y) = -l4\^\ln( y -±^] + uJ^-\-U 2 



\l + zj \ z J \\ + zJ \l + z 



H(l-z,0;y) = -lnyln I 1 - 7-^— ) - Li 2 



\l-z 



H(l-z,l;y) = ~ ln 2 (l-y) + ln(l-y-z)ln(l-y)- Li 2 (^-j + Li 2 (z) , 
H(z,0;y) = lnt/ln + U 2 (-|) , 



\ z J V z/ 

*(»,!;„) = -ln(l + ,)ln(^)-Li 2 ( T A_) + Li 2 (|±i 

H(l-z,z;y) = -\n(l - y - z)\n{^-^j + U 2 (z) - U 2 {y + z) , 

H(z,l-z;y) = ln(l - z)ln(^-^j - Li 2 (z) + Li 2 (y + z) , 

H(l-z,l-z;y) = ± In 2 f 1 - /y 



2! V !- z 

= ^ln 2 (^) • (A.15) 

The z-independent HPL i? (0, 1; y) and if (1, 0; y) are also part of the full set of 2dHPL at weight 2. 

At weight 3, a total of 64 2dHPL exist. Eight of these involve only (0, 1) in the vector of parameters 
and are therefore independent of z. The remaining 56 can all be expressed as combinations of generalised 
polylogarithms up to weight 3. Solving the product and integration-by-parts identities, the 2dHPL at weight 
3 can be expressed by a minimal set of 20 functions. This minimal set contains the two z-independent 
functions ff(0, 0, l;y) and ff(0, 1, l;y) and 18 genuinely two-dimensional functions: 

ff(0,0,l-z;y) = Li 3 



1-z 

H(0,0,z;y) = -Lis(-^) , 



F(0,l,l-,;y) = (t^) + ^3 (^^j - Li 3 ( T -^- ) + Li 3 ( y) 

- Li3 (t^) - ln i 1 - ih) U ^ - ln i 1 - ih ) Li2 ^ 

-iln 3 (l-y-z) + iln 3 (l-z) , 
6 6 

= Li3 (tTf) Li3 {jTz) - Li3 (frl) + Li3 {iTz) Li3(y) 

+ In («±f) Li 2 + In (»±£) Li 2(y) - lnzln(l + z) In ' V + * 

+ i ln(l + z) ln 2 (y + z) - ^ ln(l + z) In 2 z , 
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H(0,l-z,l;y) = S li2 (y) - Li 3 ( - ^ T ) + Li 3 I ) + Li 3 



(l-y)(l-z)J *\l-zj *\l-y 
- Li 3 (y) - Li 3 (z) - ln(l - y) Li 2 (j^j + H± - V) Li 2 (l/) + In(l - y) Li 2 (z) 

+ iln(l-z)ln 2 (l-y) , 



H(0,l-z,l-z;y) = S K 



y 



l - z 



H(0,l-z,z;y) = Li 3 ( (i - — ) - Li 3 (^-) - Li 3 (-L-) - Lfefa, + z) + Li 3 (.) 

4 hi l'^) Li 2 (z) + In f Hi!) Li 2 f-O - 1 ln(l - z) In 2 ^ + * 



1- zj 2 v ' \ z 

y( 1 + z )\ T ■ i y\ T • f 1 - y\ , T • 1 

. Li 3 



if(0,z,l;y) = -S h2 (y)+U 3 [-^^)-U 3 (-^)-U 3 [- 

+ Li 3 (y) + ln(l - y) Li 2 (-|) + ln(l - y) Li 2 ~ ln(1 ~ y) Li 2 (y) 



+ ^n(^)ln 2 (l- y )-iln 3 (l-j/) , 
H(0,z,l-z;y) = Sl , a ( (1 - ^ — ) - Sl , 2 (^-) - Sl , 2 ) + S 1|2 (y + z) 

- Si l2 (*) - Li 3 f- ^— — ) + Li 3 f-^-) + Li 3 ^ " 



(l-z)(y + z)J \y + z 

:- h, + (Li 2 („ + z)- Li a (*))-ln(»±^ Li 2 



2 / ' \ 2 7 V 1 — z 

ln(l-z)ln 2 (U±l\ + Un 2 {l-z)\J V + Z 



H(0,z,z;y) = Si, 2 (-| ) , 



i - z y v i -z 



Li,, ( I - -^1 + I.:;; (\ [ " " 



z 



H(l,l-z,z;y) = 2Si, 2 (y) - 2Li 3 — — + Li 3 ■ - - Li : 



+ Li 3 ( ^-py).) - Li 3 (-^—^-) - U 3 (y + z)+ Li 3 + 2Li 3 (z) 



y + z y V y + z 

2^ „ ,\ 1 i„.,„2/ / i-y 



+21n(l-y)Li 2 (y) + lnyln 2 (l-y)--lnzln ( , 
- In z ln(l + z) In (^±i) - 1 ln(l - y) ln 2 (y + z) + 1 In (|±£) ln 2 (l + z) 
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f z 


1 - 


z 1 


1 - 


■y 


1 + 





+ iln 3 (y + z)- iln 3 z , 
6 6 

H(l,z,l-z;y) = S x , a ( ^f^ ) ~ Si,2 (y~) ~ Sl > 2 ^ + z ) ~ Si, 2 (2 2 ) + 2 S^) 

Li 3 f TT-7 ) - Li 3 f ^) + ln (! + z ) Li 2 (t4^) ~ Hl + Z ) Li2 & 

m ( ^_ ) n 2 ( i-^) - m (i - 2 2 ) li 2 (_L.) + m Q±f ) u 2 ( y + z) 

In I ^— - I Li 2 (z)+lnzln(l-y)ln(l-z) + -ln 2 (l-y)ln(^^ 
Vi-y/ 2 \V + z 

+ ln(l - y) ln(l + z) In - 1 ln(l - y) ln 2 (l + z) 



2 V! + z / V 2 7 \ z J \l + zj \1 + Z 

y + z 



Li 3 



1 + 2 



-ff (1 — 2, 1, 1; y) = — ln 2 zln(l — z) In 2 z ln(l — y — z) + In z Li 2 [ 1 — 



2 ' 2 v " ' " V 2 

m*Li 2 I - Si, 2 fl-— + Si, 2 fl ' 



2 / V 2 / \ 2 



ff(l - 2, z, z; y) = -\ ln(l - y - 2) In 2 f ^±i) - In ( ) Li 2 (y + 2) - Lis(z) + Li 3 (y + 2) , 



2 V z / V z 

H(z,l,l;y) = iln 2 (l-y)lnf^) + ln(l - y) Li 2 (\-^) + U. 1 ' 



2 ' Vl + 2/ ' \\ + z) \l + 2 

i-y 



- Li 3 



1 + 2 



ff(«,l-a,l-*;y) = ^My + ^ln 2 ^ 1 + ln( 1 Li a (l - y - z) + Li 3 (l - z) 

-Li 3 (l-y-2). (A.16) 

To illustrate the relation of arbitrary 2dHPL at weight 3 to the respective minimal set, let us consider 
H(l — 2, 1,0; y), which is in fact the only 2dHPL appearing in the divergent 1/e-part of the planar master 
integrals. By writing out the integration-by-parts identity (A. 9) for H(0, 1, 1 — 2; y), we immediately arrive 
at: 

H(l - 2, 1, 0; y) = H(l - 2; y)H (1, 0; y) - H(l, 1 - 2; y)H (0; y) + H(0, 1, 1 - 2; y) . (A.17) 

The full basis of 2dHPL at weight 4 contains 256 functions, of which 16 are independent of 2. The 
minimal set at this weight consists of 60 functions, 3 of them independent of 2. Only particular combinations 
of 2dHPL at weight 4 can be expressed in terms of generalised polylogarithms. In general, the 2dHPL at 
weight 4 can always be written as onc-dimensional integrals over the 2dHPL of weight 3 given above: 

H(n, m 3 ; y) = f dy'f(n, y')H(m 3 , y') . (A.18) 
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A. 2.1 Limiting cases 

If y or z are fixed to particular numerical values, or if they arc related to each other, 2dHPL reduce to 
ordinary HPL. We have to derive reduction formulae for several cases, which are relevant to match the 
boundary conditions of the differential equations we study in this paper. 

In z = and 2 = 1, the reduction from 2dHPL to HPL is a mere substitution on the components of 
rh w , combined with an adjustment of the overall sign of the HPL under consideration. The corresponding 
rules can be derived by writing out the 2dHPL as a multiple integral by repeated application of (|A.5|): 



H(m w ;y) 
valid for m w ^ 0, and 



dtif(mi;ti 



di 2 /(m 2 ;i 2 ) 



dt w f (jYlw , t v 



H(m w ;y) 



dii/(mi;ti 



dt 2 /(m 2 ;i 2 



t v -i 



dt v f(m v ;t v )H(0 

W — V j tv) 



(A.19) 



(A.20) 



valid for fh w — (m v , 0^,-!,) with m v ^ 0. 

If m w z (respectively m v z in the second case), H(rh w ;y) remains finite in the limit z = 0. This 
limit is obtained by replacing /(z;i;) — > /(0;tj) and /(l — z;ti) — > in ( A.l£ , A~2(i| ) and subsequent 

integration. For m w = z (respectively m v = z), H(rh w ;y) diverges proportional to lnz in z — * 0. This 
divergence can be made explicit by the use of the HPL- algebra, such that H(m w ; y) can be written as linear 
combination of logarithmically divergent and finite terms in z — > 0. 

If m w ^ 1 — z, (respectively m v ^ 1 — z), H(m w ; y) remains finite in the limit z = 1, it develops divergent 
terms proportional to ln(l — z) otherwise. The limit z = 1 is obtained by replacing /(z; ti) —> /(— 1; ti) and 
/(l — z;ti) — > —f(0;ti) in (A.19 A.2C) and subsequent integration. 

To reduce 2dHPL in y = 1 — z, y = — z or y = 1 into HPL of argument z, one uses 



H(m(z);y(z))=H(m(z = 0);y(z = 0))+ / dz'-^H(m(z');y(z')) , 



(A.21) 



where the boundary z = can be replaced by z = 1 if H(m(z = 0); y(z = 0)) is divergent. The derivative 
■j^H(m(z);y(z)) is then carried out on the multiple integral representation of H(m(z);y(z)) (A.1S.A.2C). 
The resulting [/(m^ii)] 2 are reduced to single powers by repeated integration by parts and partial frac- 
tioning. The resulting multiple integral can then be identified as multiple integral representation of an 
ordinary HPL H(m; z). As an example, we evaluate H(l, 1 — z, y) in y = 1 — z : 



H(l,l-z;l-z) 



i?(l,Q;0)+ / dz'—H(l,l-z';l-z') 
Ji dz' 



l-z' 



dti 



1 -t 



dt 2 



dz' 



dz' 



l-z' 



1 Jo 



dU 



1 - t 2 - z' 

rl-z' 



1-h-z' 

l-z' 







dU 



d<i 

1^*1 Jo (1 - <2 



l-z' 



1 



z 



(i 



dii 



dz' 



1 j + i)F(o ;i O 



l-z' z' 
iJ(0,0;z) + ff(l,0;z) -fT(l,0;l) 



(A.22) 
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To derive the transformation formulae in y = 1 — z, y = — z or y = 1 for 2dHPL up to weight 4, we 
have programmed the underlying algorithm in FORM |21| . It is evident from the above example that the 
transformation is carried out in a recursive manner, i.e. by using the results for 2dHPL of weight w — 1 in 
the transformation of 2dHPL of weight w. 



A. 2. 2 Interchange of arguments 

To check our results for the master integrals, which were obtained in terms of H(m(z);y) by solving the y 
differential equations, we inserted them into the corresponding z differential equations. To carry out this 
check, we have to convert the 2dHPL H(m{z); y) into 2dHPL H(m(y); z) and ordinary HPL H(m; y). This 
conversion is made using 



H(m(z);y)=H(m(z = 0);y) + / dz' -^H(m(z');y) , 



(A.23) 



where the boundary z — can be replaced by z = 1 if H(rh(z = 0);y) is divergent. The z-derivative is 
again carried out on the multiple integral representation of H(m(z);y) (A. 19, A. 20). Repeated application 
of integration by parts and partial fractioning does then generate a form that can be identified as a bilinear 
combination of 2dHPL H(m(y); z) and ordinary HPL H(m;y). As an example, consider 



H(l,l~z;y) = H(l,l;y)+ f dz'-^H (1, 1 - z'; y) 

Jo dz' 

, d f v dti 



H(i,i;y) 



dz' 



dz> 



1 - 1 



dU 



dz' 



dti 



1 - 1 



i Jo 

ti 



l-h-z' 

dt 2 



i Jo 



(l-t 2 -z') 



t\2 



= H(l,l;y)+ f dz' f dh 
Jo Jo 

= H(l,l;y)+ f dz' 
Jo 



1 



1 - z 

1 

z 



1 



l-ti 
1 



z'l-U-z' 



---)H(l;y) + -H(l-z';y) 

1 — z' z J z' 

= H(l, 1; y) - H(l;y) [H(0; z) + H(l;z)} + f % [HQ. - y; z') + H(l;y) - H(l; z')\ 

Jo z 

= H(l,l;y) + H(0,l-y;z)-H{0,l;z)-H(l;y)H(l;z) . (A.24) 

Interchange formulae have been derived for 2dHPL up to weight 4 by programming the underlying algorithm 
in FORM Like the transformation algorithm presented above, the interchange algorithm also works 
recursively by using interchange formulae obtained for lower weights. 



A. 2. 3 Analytic continuation 

The 2dHPL introduced in this appendix are real for the following range of arguments: < y < 1, < 
z < 1 — y. This region of y — S13/S123 and z — S23/S123 is the kinematical region corresponding to a 
P123 — *• Pi +P2 + P3 particle decay process, as in e + e~ — > 3 Jets. To use the results obtained in this paper 
in the computation of virtual corrections to other processes, which are related to 1 — > 3 particle decay by 
crossing, one needs the analytic continuation of 2dHPL into regions of y and z outside the above-mentioned 
area. 

This task does a priori not seem feasible, since the causal prescription 

S12 = {pi+P2) 2 +iS 
s%3 = (pi+Pa) 2 +iS 

S23 = iP2+p 3 ) 2 +iS (A.25) 
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z 




Figure 1 : Kinematic regions for the analytic continuation of the 2dHPL (example) 



is broken by the elimination of S\2 = S123 — s 13 — S23 = 5123(1 ~ U ~ z ) from all expressions. The 2dHPL 
formalism does however allow to recover the correct analytic structure of the master integrals. 

In this subsection, we discuss the analytic continuation of 2dHPL to the kinematical situation relevant 
for vector boson plus jet production in hadron collisions: 2 — > 2 scattering with one final state momentum 
off shell. The kinematical region for this process, with p\ and p^ being the incoming momenta, is defined 

by 

S13 > S123 , > S 2 3 > S123 - S13 , (A.26) 

corresponding to 

y>l,0>z>l-y. (A.27) 

This kinematical region is denoted by (IV) in Figure [l], while the region < y < 1, < z < 1 — y is denoted 
by (I). To perform the analytic continuation of the 2dHPL from region (I) to region (IV), we proceed as 
follows: 

1. Continuation from region (I) to the boundary of regions (II) and (III), defined byy=l,0<z<l. 
After separation of th e diver gent ln(l — 7/)-terms by using ( A. 7 ), this continuation proceeds recursively 
as outlined in Section A. 2.1 . By specifying the continuation of 

H(l-z;y= 1) = -ln(l - 1 - z + iS) +ln(l - z) = -lnz-wr + ln(l -z) = -H(0; z) - H(l; z) -wr , 

(A.28) 

all imaginary parts of H(m(z); 1) are specified. 

2. Transformation of variables in region (III): y = 1/(1 — u) and z — v/(l — u). H(m(z);y) are then 
expressed as linear combination of 2dHPL H(m(v); u) and HPL 7J(m; v) by 

H(m{z);y)=H(m(z = v);y = l) + J du'-^H (fh (z = j^j) 5 V = Y~^) > ( A - 29 ) 
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which is well defined since all divergent ln(l — j/)-terms were separated off beforehand. These terms 
are continued according to 

H(l;y) = -]n(l-y-iS) = - ln(-u-iS)+ln(l-u) = - lniH-wr+m(l-ii) = -H({):u)-H(l:u)+in . 

(A.30) 

In region (III), one has < u < 1 and < v < 1 — u, such that the H(m(v); u) and H(rh; v) are real. 

3. Interchange of arguments: The H(m(v);u) are transformed into a linear combination of H (m(u); v) 
and H(m; u). All terms proportional to lnu are separated off using the product identities, such that 
the remaining terms are finite in the limit v — 0, which is the boundary to region (IV). 

4. Transformation of variables in region (IV): u = 1 — r and v = — £. H(m(u); v) are then expressed as 
linear combination of 2dHPL i/(™(r);^) and HPL i?(m;r) by 



H(m(u);v) = H(rh(u);v = 0) 



l ^'^ H ( ™ {U = 1 ~ T) ;V = ' (A - 31) 



which is again well defined. The remaining lnw-terms are continued according to: 

H(0;v) = In [ ] = ln(-£ + W) = ln^ + iTr = ff(0;£) + in . (A.32) 



y 

The H(rh; u = 1 — r) are transformed into H(rh; r) by using the standard HPL formulae of 
The newly introduced variables fulfil 

< T = £Hi < i ) an d o < £ = — ^ < 1 - r , (A.33) 

Sl3 Sl3 

which is the region where the 2dHPL H (to(t); £) and the HPL _ff (m; r) are real. The imaginary parts have 
been separated off, and appear explicitly. 

With this example, we have demonstrated how the analytic continuation of 2dHPL can be carried out 



in practice. The main point to notice is that the analytic continuation formulae (A. 28), (A.30) and (A.32) 
appear to be mutually inconsistent as far as the assignment of imaginary parts to the variables y and z is 
concerned. This apparent inconsistency is however only an artifact of the identification s\2 = si23(^ — y — z). 
Considering the HPL and 2dHPL at weight 1 as the basic objects for analytic continuation (and not the 
variables y and z) this inconsistency is no longer present. It is sufficient to specify the analytic continuation 
at weight 1, the continuation for higher weights is obtained recursively from the lower weights. 
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